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Abstract 

The supersymmetric extension of the standard model with an ad- 
ditional gauge singlet is analysed in detail in the light of the recent 
experimental bounds on supersymmetric particles. The useful part of 
the parameter space and the particle spectrum are displayed. We find 
that once the recent bounds on the chargino mass are imposed, all 
other new particles practically satisfy the present experimental limits. 
Special attention is given to particles to be searched for in the future 
experiments. The singlet fields tend to decouple and give rise to an 
effective MSSM, enlarging the validity of many phenomenological anal- 
yses based in the minimal field content. However, in some ranges of 
the parameters the singlino is the lightest neutralino, which modifies 
the signature for susy particles. Simple analytical approximations are 
developed that qualitatively explain the numerical results. 
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1 Introduction 



The addition of a SU(3) x SU(2) x U(l) gauge singlet supermultiplet to 
the MSSM (minimal supersymmetric extension of the Standard Model) has 
been motivated by several issues: (i) the solution of the so-called /i-problem; 
(ii) the possibility of spontaneous breaking of the CP symmetry; (hi) the 
upper limits on the Higgs mass that are relatively low in the MSSM could 
be relaxed; (iv) the neutral particle spectrum is enriched by two scalars and 
one fermion, which could mix with other Higgs bosons and neutralinos and 
modify their physical properties; (v) possible deviations from the MSSM 
could urge to enlarge the experimental searches, which have been centred 
upon the MSSM. We refer here to the supersymmetric model with a gauge 
singlet in its simplest version, as defined below, as the (M+1)SSM. 

However, according to our analysis, the strong constraints from LEP ex- 
periments, when combined with theoretical and phenomenological require- 
ments, select the (M+1)SSM parameters in such a way that the gauge singlet 
scalars and fermions are basically decoupled and the rest of the spectrum 
effectively reproduces a restricted class of the MSSM. Because of this sim- 
ilarity, the motivations (ii)-(v) fade out as these issues are excluded or di- 
minished after such a detailed investigation. Nevertheless, this also means 
that phenomenological analyses of the MSSM apply to a large extent to the 
(M+1)SSM as well! Still, if supersymmetry is found, one can go further and 
look for some peculiarites that could distinguish the two kind of models, 
specially in the case in which the lightest supersymmetric particle is the sin- 
glet fermion, or singlino. The (M+1)SSM can also be experimentally tested 
from the fact that the effective parameters in the equivalent MSSM are very 
constrained. 

The present paper is devoted to a thorough discussion of the model in its 
version with only dimensionless couplings in its superpotential, and univer- 
sal boundary conditions for the soft supersymmetric breaking terms at the 
GUT scale. With these assumptions, the dimensionality of the parameter 
space of the (M+1)SSM equals that of the MSSM. Our aim is to under- 
stand the pattern of the particle spectrum, correlations among the particle 
masses, and features which are in common with or different from the MSSM. 
An analysis of the complete parameters space of the model requires numer- 
ical methods, as a scanning over the independent parameters at the GUT 
scale, their subsequent evolution using the RGE down to the weak scale, the 
minimization of the Higgs potential, the implementation of theoretical and 
experimental constraints and finally the computation of the particle masses 
and some mixing angles. 

Many phenomenological analyses of the model have appeared before 
0, ||, H|, ||, ||]. The purpose of the present work is twofold. First, we imple- 
ment the most recent experimental constraints on sparticles, Higgs and top 
quark masses from the LEP and the Tevatron [||, [7|]. Second, we explain the 
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results of our numerical procedure with the help of analytic approximations 
and inequalities among the parameters. Once the vacua which break colour, 
electromagnetism or correspond to tan/3 = oo are forbidden and the ex- 
perimental constraints are imposed, these formulae allow us to understand 
the essential features of the bounds on and correlations among some par- 
ticle masses. Quite a lot of attention is paid to the property found in 
namely, that the singlet scalar acquires a v.e.v. which is much larger than 
the Higgs ones. Indeed, this has the important consequence that the singlet 
fermion and scalars quite decouple from the the other multiplets as already 
mentioned. 

In section |2| we review the motivations for the model, its essential differ- 
ences with respect to the MSSM due to the enlargement of Higgs sector and 
the status of the spontaneous CP violation. The universality assumption 
for the soft terms and some of its consequences are explained in section ||. 
Section [|] contains the minimum conditions for the scalar potential as well as 
the scalar Higgs masses with the radiative corrections included in an approx- 
imation that is generally good, though we use the full expressions || in our 
numerical calculations. In section [| we derive some approximate constraints 
from the stability of the vacuum to understand the pattern of the physical 
predictions. In section || we discuss some correlations in the mass spectrum 
along the same approximations. Then we summarize our implementation of 
the experimental and theoretical constraints on the model, in section 0. The 
main results of the careful scanning of the parameter space are displayed in 
the figures in section |8|, with some emphasis on the sparticles more likely to 
be discovered at the colliders, with the exception of the Higgs sector that has 
been already presented (5[. The results are summarized in the last section. 
Some useful solutions of the RGE for the soft terms (with generic boundary 
conditions, for completeness) are listed in the Appendix. 

2 General properties of the (M+1)SSM 

The superpotential of the (M+1)SSM (where family mixing is neglected and 
irrelevant for the analysis in this paper) is given by 

W = XSH^ + ^S 3 + h t TH 2 Q 3 + h b BH 1 Q 3 + h T EH 1 L 

+ (similar terms for the lighter quarks and leptons) (1) 

where the chiral superfields are denoted as follows: S (singlet), H\ and H 2 
(Higgs doublets), T (antitop), B (antibottom), E (antitau), Q 3 (top-bottom 
doublet), L (tau- neutrino doublet). Let us briefly review the present status 
of these motivations for the (M+1)SSM. 
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A supersymmetric mass term fiHi ■ H2 is added to the MSSM superpo- 
tential in order to avoid a potentially dangerous Peccei-Quinn symmetry and 
to allow for an acceptable minimum of the Higgs potential. The analysis of 
the SU(2) x U(l) breaking then requires // to be roughly of the same order 
of magnitude as the soft terms that embody the supersymmetry breaking 
effects in the two Higgs sectors. This interesting hierarchy problem, the 
so-called //-puzzle, has motivated many a suggestion to link the supersym- 
metric mass // to supersymmetry breaking H. In any instance, some kind 
of new physics is requested and in this sense the MSSM that looks mostly 
economical at low energies is not necessarily so at the unification (GUT) 
level. A natural solution [1C] is provided by the addition of a singlet S. 
Soft terms that break supersymmetry can induce a v.e.v. for S and an ef- 
fective supersymmetric mass [i = XS of the right order magnitude. This 
mechanism asks for the kS 3 /3 self-coupling in the superpotential which also 
breaks the unwanted Peccei-Quinn symmetry. However, this coupling would 
be forbidden and the singlet model inoperative if S has non-trivial quantum 
numbers with respect to broken symmetries beyond those of the Standard 
Model. This phenomenon is frequently found in superstring compactifica- 
tions. We take here the pragmatic viewpoint that all the dimension three 
interactions that preserve the low-energy local and global symmetries are to 
be included. As we shall discuss later, very low values of the coupling k are 
consistent with the phenomenological constraints. In this case, it could be 
a relic of non-renormalizable interactions after the decoupling of, e.g., extra 
U (1) gauge symmetries. Actually, the "sliding" singlet v.e.v. compensates 
for the smallness of the couplings n and A and this considerably reduces the 
fine-tuning problem of the //-term. Such a possibility is indeed supported 
by our analysis. 

Supersymmetric mass terms are forbidden by the Z% symmetry of the 
purely cubic superpotential assumed in (||). This discrete symmetry also 
prevents S from getting a large v.e.v. through large radiative corrections to 
singlet tadpole terms, which invalidate other supersymmetric models with 
gauge singlets |TT| , |l2|| . On the other hand, the ^3 symmetry can give rise to 
cosmological problems in the form of domain walls which are produced dur- 
ing the electroweak phase transition [13]. Possible solutions to this problem 
have been discussed in [13| and include inflation of the weak scale, embed- 
ding of the discrete symmetry into a gauge symmetry at the Planck scale 
or re-introducing of the //-term in the superpotential. (In the presence of 
the singlet, // can be orders of magnitudes smaller than the weak scale, in 
which case our subsequent analysis is not affected; however, the //-problem 
was part of the motivations of the model.) In any case, since we concentrate 
in this paper on the particle physics aspects of the model and do not dis- 
cuss cosmological issues, our results will be independent of the form of the 
solutions of the domain wall problem. 

The study of spontaneous CP breaking in the (M+1)SSM also differs 
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from the MSSM. Indeed, CP violating vacua can be induced only if there 
are more than one phase-dependent interaction in the scalar potential. In 
the MSSM, the only tree-level term is the soft scalar mass term B^iHiH^- 
Analytic terms of higher dimensions generated at the one loop level can con- 
spirate with the tree-level one to induce spontaneous CP violation. However, 
this has been shown [14] to entail a very low upper limit on one of the scalar 
masses which is inconsistent with present experiments at LEP. The classical 
scalar potential processes a richer structure in the (M+1)SSM, with three 
additional complex scalar couplings (one of dimension four, two cubic soft 
terms, see next section). The conditions for CP spontaneous breaking in 
this type of models have been shown [15] to imply a negative eigenvalue 
in the scalar (mass) 2 matrix. However, top quark-squark loop corrections 
could shift this eigenvalue to positive values. Thus spontaneous CP viola- 
tion would occur for appropriate values of the parameters and one or two 
light scalars are predicted |l(|. This issue certainly deserves a more detailed 
comparison with phenomenology; it lies outside the purview of the present 
investigation. Because the very small region in the (M+1)SSM parameter 
space consistent with spontaneous CP violation and those selected by our 
CP conserving solutions are apart, we disregard ah initio the question of 
non-trivial phase vacua. 

Upper bounds have been derived for the mass of the lightest Higgs scalar 
in the MSSM that take into account the large one loop corrections [[l?]] as 
well as some two- loop effects [jDf| . This prediction is of paramount impor- 
tance in view of the experimental searches of the Higgs scalar in the LEP 200 
and the LHC. The (M+1)SSM has an additional contribution cx |A| 2 to the 



Higgs (mass) 2 . It has been shown M, 0, 19] that requiring the theory to re- 



main perturbative up to the GUT scale reduces this additional contribution 
to < O(30 GeV). The upper bounds on the lightest Higgs boson mass are 
considerably lowered if one also imposes all phenomenological constraints 
needed for the model to be realistic. This has been done in the case of uni- 
versal supergravity couplings [||, ||] and also for some special non-universal 
models |20[| . Note that a relatively small increase in the Higgs mass may 
affect the experimental signatures. Since the predictions of the phenomeno- 
logically constrained (M+1)SSM for the Higgs sector has been the subject of 
a recent publication || , we shall only give a short account of this important 
matter here. 

As for the particle spectrum, this model presents several new features 
as compared to the MSSM. The supermultiplet S consists of two scalars of 
opposite CP and one Majorana fermion. After the SU{2) x U(l) breaking, 
these states get mixed with the neutral states in the Higgs doublet and gauge 
supermultiplets, leading to a slightly richer spectrum. In particular, since 
the visibility of physical states depends on their couplings to the gauge 
bosons and to matter, some particles could be relatively light if they are 
mainly gauge singlets. Thus, some of the limits on neutral particles obtained 
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at LEP in the MSSM framework do not immediately apply to the alternative 
(M+1)SSM. The discussion in this paper as well as in our previous ones [§, ||] 
is also aimed to persuade experimentalists to analyze their data on a more 
general basis. 

In the k — > 0, A — > limit, the gauge singlet decouples and the remaining 
theory looks like the MSSM However, this limit leads to a particular 
version of the MSSM with a strong correlation among the soft parameters, 
as we shall discuss later on. The singlet scalar takes a very large v.e.v. of 
0(1713/2/ k), but the fermion singlet and both scalar singlets are rather light. 

Of course, the large dimension of the scalar field manifold allows for a rich 
structure in the overall scalar potential and leads to the existence of many 
local extrema pl| . One has to carefully check that the phenomenologically 
viable vacuum is favoured with respect to the unwanted ones. Among these 
there are the usual SU(3) C x U(l) cm breaking ones, corresponding to vacua 



with quark and lepton quantum numbers |p2| , 21|| , [23] . In the (M+1)SSM 
there is the additional possibility of a non-vanishing charged Higgs field vac- 
uum. However, the necessary condition for this particular charge violating 
classical solution is a relatively large coupling of the Higgses to the singlet 
[24], so that this situation is quite naturally avoided by a bound on that 
coupling. Our theoretical and experimental constraints turn out to auto- 
matically selected the right region of the parameter space. 

Actually, as we shall discuss here below, an acceptable minimum of the 
potential requires some inequalities among the soft supersymmetry breaking 
terms which would tend to induce SU(3) C x U(l) cia breaking vacua as well. 
The radiative corrections due to the gauge interactions then allow for a com- 
promise while restricting the useful region in the parameter space. Since the 
electroweak gauge symmetry breaking is induced by radiative corrections in- 
volving the top Yukawa coupling - a well-known property in supersymmetric 
extensions of the Standard Model - the existence of a physically acceptable 
vacuum in the (M+1)SSM is a purely quantum effect. 



3 Soft terms and boundary conditions 

The supersymmetric part of the scalar potential is obtained from the super- 
potential through the well-known expression as a sum of F and D terms. 
The supersymmetry breaking part of the lagrangian contains all soft terms 
consistent with the symmetries of the superpotential: gaugino masses, trilin- 
ear analytic scalar couplings and scalar mass terms. In the (M+1)SSM they 
are as follows (with the same notation for the chiral matter supermultiplets 
and their first component complex scalars) : 

Asoft = (M1A1A1 + M 2 A 2 A 2 + M3A3A3+ 
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+A X XSH 1 ■ H 2 + A K ^S 3 + A t h t Q 3 ■ H 2 T C + h.cj + 

+m\ l^l 2 + m\ \H 2 \ 2 + m\ \S\ 2 + m 2 Qa \Q 3 \ 2 + mf, \T C \ 2 + 
+... " (2) 

where Ai, A2, A3 are the gauginos associated to the U(1)y, SU(2) and 577(3) 
gauge symmetries. Only the terms involving the top scalars have been re- 
tained in (||), with Q3 = (T, B), but similar interactions for all other squarks 
and sleptons are to be understood. All the parameters are taken to be real 
(up to irrelevant phases) and flavour mixing is neglected. Notice that these 
soft terms contain all the dimensionful parameters of the model, which are 
all naturally related to the supersymmetry breaking scale M susy . In order 
to study the particle spectrum and 577(2) x 17(1) symmetry breaking in the 
(M+1)SSM, one should take the parameters in (|l|) and @ at the relevant 
scales, the Fermi scale v = 174 GeV or M susy . Instead, it is natural to fix 
these parameters at the level of the underlying supergravity theory. We 
shall consider here the unification scale Agut ~ 10 16 GeV to define them. 
This is just a few orders of magnitude below the Planck mass and one can 
presumable assume an effective renormalizable theory at Agut- Once the 
parameters are given at AquT; their corresponding values at v or M susy are 
calculated by the integration of the renormalization group equations (RGE) 
for the (M+1)SSM plf . We present in the appendix the analytic solutions 
of the (M+1)SSM RGE in the approximation h 2 S> A 2 , k 2 , which turns out 
to be good for most of the set of parameters that will be selected on a 
phenomenological basis. However, for the sake of precision, we numerically 
integrate the RGE in our analysis, the expressions in the Appendix being 
reserved for analytical discussions of the numerical results. In practice, we 
renormalize the parameters down to the scale v through the supersymmet- 
ric RGE, and then apply the Coleman- Weinberg corrections to the Higgs 
potential to take into account the decoupling of the stops at M susy . 

The important hypothesis to be made in this paper is the so-called "uni- 
versality" or flavour-independence for the soft-terms in (||). For the super- 
symmetry breaking scalar interactions it amounts to assume a common value 
for all the parameters Ai in (||) (at Agut); 

A x = A K = A t = A b = A T = ... = A , (3) 
and the same values (at Agut) for the scalar masses in (|J), 

2 2 2 2 2 2 f A\ 

m 1 = m 2 = m s = mg 3 = m T = ... = m . (4) 

This universality follows if the direction of supersymmetry breaking, char- 
acterized by the goldstino components, corresponds to fields with purely 
gravitational couplings to the relevant chiral superfields (for instance, if it 
is along the dilaton-axion direction in superstring inspired theories [p6[ ). Of 
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course, it follows in general if the goldstino components have equal couplings 
to all relevant chiral superfields; an example is given by the "large radius lim- 
its" of orbifold compactifications in string theory if supersymmetry breaking 
occur along the modulus direction |ffif| . 

As for gaugino masses, the universality hypothesis at Agut reads 

Ah (A GUT ) = M 2 (A GU t) = M 3 (A GU t) = M . (5) 

In grand unified gauge theories many of these relations will be required by 
the enlarged gauge group at the tree-level. At the one loop level one expects 
some threshold corrections to the universality of the soft terms already at 
the reference scale Agut- The deviations from (||) - (0) in the matter sector 
can be phenomenologically important because of the existent experimental 
bounds on FCNC effects but they are less relevant to our discussion in this 
paper. 

Deviations from the universality assumption are also possible in the 
framework of superstring theory, since the moduli sector can couple dif- 
ferently to the various chiral fields. Since the assumptions ^4 k (Agut) = 
A T (A GVT ) and mf(A G uT) = m|(A G uT) = m|(A G uT) = m|(A G uT) = 
m|(AGijT) play & n important role in our phenomenological analysis, some 
of our results depend on the universality of the supersymmetry breaking pa- 
rameters. However a similar analysis of the (M+1)SSM has been published 
[20 1 for the case of non- universal soft terms in the framework of an orbifold 



compactification of the superstring, and these results are quite analogous to 
those presented here. 

The most important consequences of the flavour independence of soft 
terms concern the SU(2) x U(l) symmetry breaking. It could be induced 
by a negative value of m|, but for the potential to be bounded from below 
and to avoid colour or e.m. symmetry breaking the universal parameter mp 
has to be positive, so that toKAqut) > 0. For this reason, < has 

to be enforced by renormalization effects. This is a well-known feature of 
the MSSM as well, entailing a lower limit on the top Yukawa coupling. In 
the (M+1)SSM, the sign of m\ could be made to flip by the coupling of 
the singlet fields even for a light top quark [21]. The experimental value of 



top quark mass is perfectly consistent, however, with a negative m| at low 
energies from the effects of top-stop quantum loops. 

Furthermore the cubic soft terms can also induce spontaneous symme- 
try breaking if the value of the parameter A in (||) is large enough. Hence 
large values of A K are preferred to give the singlet scalar a v.e.v.. However, 
with the universality assumption (|3|) , the squark and sleptons could also get 



non-vanishing v.e.v. 's [22, 21, 23|. Interestingly enough, the gauge interac- 
tions reduce the ^4-parameters for matter fields while preserving the value 
of A K , as can be seen from the expressions in the Appendix. Thus, phe- 
nomenologically consistent (M+l)SSM's with the condition (§) , @ , © , can 
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be build at the price of restricting the space of the supersymmetry breaking 
parameters. This is the most important restriction imposed on the model 
by the universality hypothesis, as compared to the MSSM J2|. Otherwise, 
the effects of this universality on the SU(2) x U(l) spontaneous breaking 
are similar to those in the MSSM. 

Another consequence of this assumption is to reduce the number of free 
bare parameters in the (M+1)SSM to five (before imposing the physical 
constraints) just like in the MSSM: the dimensionful parameters B and p of 
the latter are replaced in the former by the dimensionless parameters A and 

K. 

At low energies, the soft terms as obtained from the RGE depend on the 
parameters Mq, Aq, m 2 ,, on the choice of the GUT scale, Agut> through t = 
ln(AGUTA') / (l67r 2 ) (a typical value being t ~ .21), and on the gauge and 
Yukawa couplings. This calculation has to be done numerically in general. 
However, from the dimensionality of these running parameters, they can be 
generically expressed at low energy in terms of their universal initial values 
as 



Ai(t) = aiAo + biMo 

m?(i) = Ziml + XiAl + yiAoMo + WiM^ (6) 

where the coefficients depend on the various gauge and Yukawa couplings. 
For the simple case where the top Yukawa coupling dominates over the 
others, the dependence on the latter reduces to simple expressions in terms 
of the ratio p = h 2 (v) / 7i 2 rit where h cvlt ~ 1.1 is the fixed point value of the 
top Yukawa coupling. The results of this approximation are presented in the 
Appendix and are used for the sake of some qualitative understanding of the 
numerical results. With the recent measurements of the top quark mass, p 
is bounded by p > 2/3 (this is higher than the minimum value required to 
induce mediative SU(2) x U(l) breaking). As an example, we list the low 
energy soft parameters at the fixed point value p = 1. The coefficients are 
approximated to a reasonable precision. 

A K ~ A A x ~ - (A) - M ) A t ~2M 

m i — m o + m 2 — ~~Y ~ 3Mq m 2 s ~ Trig 

m\ ~ 4M 2 m 2 Q3 ~ -m 2 , + 6M 2 (7) 
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4 Mass spectrum in the Higgs sector and vacuum 
structure 

Let us first restrict the discussion to the physically relevant vacua with (real) 
v.e.v.'s for the field (H®) = h,2,(H®} = hi,(S) = s. The scalar potential 
restricted to this sector reads (g 2 = g\ + <?|) 

V = { K S 2 + \H°H$) 2 + \ 2 \S\ 2 (\H° 1 \ 2 + \H°\ 2 ) + 

+ (a x \SH°H$ + Ak^S 3 + /i.e.) + mi \H^\ 2 + m\ |tf 2 °| 2 + 

+™ll^l 2 + ^(|F 1 °| 2 -|^| 2 ) 2 + y rad (8) 

where V Ta d is the quantum correction to the effective scalar potential beyond 
the RGE effects included in the running parameters of (||). These radiative 
corrections are important, in particular in the evaluation of the lightest 
Higgs boson mass. The relevant contribution comes from the top-stop sector 
because of the relatively large value of the top Yukawa coupling. Hence we 
only include the top quark and the two stop states, T\ and T 2 in the usual 
Coleman- Weinberg expression, 

1 TTX^ 

V ™ d= 64^ 5Trm4ln Q2 ( 9 ) 

The appropriate mass terms in (^) are then the top quark mass m 2 and 
the stop masses 

m Tl 2 = m t + — (rriQz + m T ) ± W, where 

W 2 = j[rn 2 Q3 +rri^) +h 2 \A t h 2 + Xsh^ 2 , 
m t = h t h 2 . (10) 

All the parameters in (^) and @ are to be taken at the scale Q 2 ~ 
O (h 2 + h 2 ). 

These radiative corrections are mostly due to the splitting between the 
quarks and the squarks. However, even if one starts from universal pa- 
rameters at AquTj the RGE running introduce some asymmetry in the stop 
sector as can be seen from the formulae in the Appendix. The resulting 
splitting in m\, — produces a sizeable effect in V^ad an d in the scalar 
mass spectrum. The full one loop expressions for the mass matrices and the 
potential minimization are available in the literature. All the numerical cal- 
culations in the present paper are performed by using the formulae in ||. 
However, for the sake of the analytic approximations discussed below we use 
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a simple approximation by keeping only the lowest relevant order in the de- 



rn 



T 2 ) 



velopment in terms of the stop splitting, (to^ — to 2 -^) I ( m Ti + 
minimization of the scalar potential amounts to the following conditions 



The 



where 
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+ j3 t [i 2 L 



ml + A 2 (l + f) - \t (hj - hi) 
+BAfe + Pt [{rn 2 T + m 2 Qa + A 2 ) (L + e) 

- (m^ + m Q^j + 2m. 2 L 



hih 2 I 
— 1~-r 



3h 2 



= 



1 16vr 2 ' 
and the variables 



In 
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Q 2 ' 



L = In 
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o 

mf 



(11) 



(12) 



/i 
5 



As 



^4a + f + PtM (L + e) 



(13) 



have been introduced. An one loop term has been conveniently included in 
the definition of the effective parameter B (of the MSSM). 

The 3x3 mass matrix for the CP even scalars (in the same approxima- 
tion) has the following matrix elements 



to 2 ! 



"1-22 



m 2 ss 



fhi -Bfi[^j+ 4f3 t m 2 L + t A t (^2m t Z - A t ^- 

hih 2 P 2 ( h\ \ a z2 
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m 



2S 



2jr[ — ) + [B + v)ii[ — )+ fi[ — )p t [ mtZ - A t — 



z- 



m 



12 



-fh x h 2 + 2p 2 ( ^] + Bp + (ify (m t Z - At%- 



m t (A t + p cot (3) 



(14) 



The parameter L characterizes the top-stop mass splitting and can be rel- 
atively large, while e depends on the low energy scale chosen to define 
the running parameters — we take Q 2 ~ [h\ + h^j in this paper so that 
|e| ~ I In Ttf I < 1. 

The CP odd scalar (mass) 2 matrix with the would-be Goldstone boson 
projected out is as follows, 



/ a . h?+h 
- B V-fcjT 2 

V 



B-3v) fi 



(b - 3.) A^P 



3A K u- (B + 3v) fi 



. h\h,2 



\ 



(15) 



The minimization of the potential fixes Bfi and A K v to be negative, so that 
the diagonal entries in (^) are positive. The (mass) 2 of the charged scalar 
reads 



m 



1 



A 2 
9\ 



+ 



Bp, 



h\ + h 2 2 



hi h 



l"-2 



(16) 



The radiative corrections are incorporated into the CP odd neutral and 
charged scalar masses through the shift of the B parameter as defined in 

(ED i- 

Although all the numerical results in this paper are obtained by nu- 
merical minimization of the complete one loop potential (|8[), (^) and the 
numerical determination of the mass eigenvalues, it is interesting to get 
some insight into the results through simple analytic approximations. In- 
deed, many features of the particle spectrum of the (M+1)SSM are related 
to an important property, first noticed in ||], namely: the only physically 
acceptable solutions of the vacuum equations ([ll|) are large singlet solutions: 
s > ITeV S> h,2i h\. Indeed, solutions with s < h 2 lead to one or more light 
states in the spectrum that are excluded by experiments, specially by LEP 
experiments. Examples can be found in Q, where the parameters have been 
chosen to yield light spectrum. The general character of this feature of the 
(M+1)SSM has been carefully checked in our numerical analysis. Having 
settled this fact one can go on to establish several approximate consequences. 
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(A) The S v.e.v is approximately obtained from the third equation in 
(0) with h 2 = hi= 0, 



S * (17) 

The potential at this value of s has also to be negative in the allowed region of 
the parameter space, as discussed in more detail in section [|. This additional 
condition requires 

9m| < Al (18) 

This approximation turns out to be very accurate. The parameters A K and 
rm? s are only slightly renormalized if k, A <C 1. Actually, these relatively 
low values of the couplings are also required by our numerical study of the 
model, so that the condition ( |lg| ) approximately applies to the bare universal 
parameters Aq^itiq, yielding the strong constraint in the parameter space: 

Al > 9m 2 (19) 

The S v.e.v. in ( |17| ) defines the effective parameters B and fi analogous 
to those often denoted B and /i in the MSSM. Indeed, up to corrections 
of O (h^/s 2 ) , the first two equations in (11) are the same as in the MSSM 
with the introduction of these parameters. Notice that ft ~ XA k /2k should 
be of O (w.3/2) • Hence the well-known "fine-tuning" on the MSSM free 
parameter, fi ~ O (^3/2) ; is replaced here by the requirement of no strong 
hierarchy between the A and k Yukawa couplings, as the SU{2) x U(l) 
invariant scale of the singlet field is dynamically fixed by supersymmetry 
breaking, accordingly to (|l7|). 

With tan/3 = /12//11 one recovers, in this approximation, the MSSM 
relation 

tan* (, = . ™j + a/^ + A 2 + MS/2 (2Q) 

ml + (m 2 Ti + m\ 2 + A 2 ) (3 t L + jj 2 + M|/2 

Only the dominant terms of the radiative corrections are retained. Be- 
cause the denominator in (|2^) has terms of either sign (m\ < 0) it is rel- 
atively sensitive to radiative corrections which tend to decrease tan (3 with 
respect to the tree-level approximation. 



(B) In the neutral scalar mass matrices (14,15) as well as in the neu- 
tralino mass matrices, the mixing of the singlet fields to the others is always 
proportional to h\/s and hil s, hence small. This is the reason why, in 
most of the allowed points of the parameter space, the singlet sector of the 
theory, after dynamically producing the effective parameters (|l3|), is almost 
decoupled from the rest of the theory that resembles to the MSSM. 
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(C) The MSSM is formally obtained in the limit s — > oo with v and fi 
fixed (k, A <C g 2 ) ■ In this limit the singlet sector decouples from the doublet 
Higgs sector, and the singlet masses are: 



~ 8mg (\A \ + \J Al~ 8m%) (CP even scalar) 

\A \ (\A \ + yjAl -8m^j (CP odd scalar) (21) 



3 
4 



Aq\ + \f Aq — 8t?i,q I /2 (singlet fermion, singlino). 



The singlet particles become light in the small Aq limit. We argue in section 
[| that this limit is always associated to the singlet decoupling limit 
(k — > 0, s 2 3> h 2 + ^1) • Some of the physical consequences of a light singlino 
are discussed in section ||. 

The neutral and charged Higgs0 sector, as well as other neutralinos and 
charginos follow the MSSM pattern. In this sense the (M+1)SSM reduces to 
the MSSM in this limit but with an important restriction in the parameter 
space needed to dynamically implement the physically acceptable vacuum. 



5 Constraints on the parameter space from vac- 
uum stability 

The spontaneous breaking of SU(2) x U(l) symmetry in the (M+1)SSM 
restricts even more the (Mo, mj, ^4o) parameter subspace than in the MSSM 
under the same assumption of universality. The minimum discussed in sec- 
tion 3 has to be cosmologically flavourless and stable with respect to other 
possible minima. The first constraint comes from the fact that the S v.e.v. 



is induced through the "A-mechanism" [26, 10] with Aq > 9mQ, as given in 



(19). It is well known that large Aq values may induce charged and coloured 



vacua. 

In order to avoid slepton v.e.v. 's one has to impose 



A 2 C < 3 (m| + m 2 L + m\) (22) 



As can be seen from ( |lq ) there is a negative contribution A M z /~g , to m H+ . This 
is the origin of the condition A 2 < ~g 2 which is sufficient to avoid cm. charge breaking 
through (H + ^ 7^ 0, which is satisfied anyway by the models that are consistent with the 
other phenomenological constraints. 
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at a scale A < (A c /h c ) , where (A e /h e ) characterizes the scale of the would- 
be slepton v.e.v.'s. From the RGE solutions in the appendix this implies 
(t e = In (h e A GVT /A e ) /16vr 2 ) , 

{A + 2M t e (3sf + 30?)) 
< 9[m 2 + 2g 2 t c (l-g 2 t c )M 2 + 2g 2 t c (l-llg 2 t c )M 2 ] . (23) 

The most dangerous v.e.v.'s comes from the electron sector with the 
smallest Yukawa and the largest scale A e /h e ~ (l0 8 GeV^), and (22) at this 
scale finally becomes 



(A) + AM Q ) 2 < 9m 2 + 2.4M 2 . 



(24) 



This can be combined with the condition Aq > 9m 2 to obtain restrictions on 



'o 



the ratio Aq/Mq. For this purpose we introduce a parameter to characterize 
the fine-tuning of m 2 , with respect to Aq, 



a 



A 2 



9m,Q 



At 



and obtain from ( f24| ) the upper bound on a 



M 



An 



> .18 (Vl4a + 1 + sign(A)/M )) 



For Aq/Mq > 0, (g4j) yields the limit 

A 2 



o 



(25) 



(26) 



(27) 



For Aq/Mq < 0, |Ao | can become much larger than Mq at the price of a fine- 
tuning of ttiq/Aq corresponding to a of O (Mq/Aq) . Though these conditions 
are only valid under the universality assumption for soft terms, analogous 
ones can be derived along the same lines for other models (see, e.g., |20|| ). 
(Strictly speaking, ( p4[ ) can only be analytically derived (2l[| if the scalars 
are relatively degenerate. But recent numerical studies suggest that it is 
approximatelly valid under more general conditions p5| ) 

The next step is to avoid minima with (Hi) = 0,i.e., with tan/3 — > 
oo. In this analysis the radiative corrections discussed in section 3 can be 
relevant. Of course, in our numerical analysis they are included, but in 
order to understand the origin of some of the constraints on the (M+1)SSM 
parameter space, they are neglected in some of the expressions here below. 
In the effective two-step minimization of the scalar potential, one finds for 
the phenomenologically acceptable solution discussed in the previous section 
the approximate value at the minimum 
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— A 4 M A 
Knin = — iC(a)--4cos 2 2/3 



C(a) = (i±4±=) (^±^i) (28) 

where the first term comes from the pure singlet sector and the last one 
from the effective MSSM potential defined with the parameters given by 
(|i~3|). This physical minimum has to be compared with another one, with 
(S) = (Hi) = and (-fff) = — m i/2ff 2 - The potential at this minimum is 

^min = =2 (^) 

In principle one should look for the cosmological formation and stability 
of the minimum ( [28] ) with respect to (p^). For simplicity we replace the 
rigorous constraints by the simpler condition that V m \ n < V^ in , i.e., that the 
absolute minimum is the physical one. This stability condition differs from 
the corresponding one within the MSSM in two respects: (i) the presence 
of the pure singlet negative term in (pq ) favours V m i n ; (ii) in the MSSM, 
(— m 2 ,) 2 is replaced by (— m| + fi 2 ) 2 in (|29|), so that V^ in is more dangerous 
in the (M+1)SSM. The comparison of (ES) and (|29|) leads to the condition 
(at the scale v), 

|^C(a) + ^cos 2 2/?>(m 2 ) 2 . (30) 
On the other hand, from the RGE solution in the appendix one gets 

2 . / : V' ,\ 2 (, 1 



-™2> ^y-lJm^+^4p--jM ^ (31) 

and the measured values of the top mass does not allow for small p = 
h 2 /h 2 rit so that p > 2/3. In addition, in the next section we derive from the 
experimental limits on the chargino masses a bound on Mz/Mq, 

2M 



^_ >M z =g^h{ + hi . (32) 

(At this point, one is introducing experimental information on the spectrum. 
As already stressed, it is the interplay between the LEP bounds on sparticle 
masses and the vacuum conditions that constrains the parameter space.) 
Now, from (||) and ©, © can be well approximated by 

_J_ (Ml) (4P-*) m 

±VSK \A 2 J y/C&) 
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Notice that the necessary condition C(a) > 0, which follows from (|32|) and 
implies a > 0, namely, the condition fli~9|). 
For Aq/Mq > one obtains from (|26|) and (|33|) the bound 



- < 7 1 n , (34) 
9 5(4p-i) 

so that k is always relatively small, n < g/10 = 0.05. This is related to the 
large values taken by the S v.e.v.. Indeed by replacing ( pTTj ) and (32) in (|33|), 
we obtain 

-^^(Sp-l)^. (35) 

Of course, the limits ([34|) and ( |35| ) are relatively loose bounds and most 
of the solutions of our numerical studies correspond to smaller (larger) values 
of k (s, respectively). Also to be noticed is the behaviour of s for fine-tuning 
of the parameter a in (p5|): s/ \Jh\ + h% — > (4p — 1/2) / y/a. as a — > 0, while 
k decreases as -y/a. 

For Aq/Mq < 0, the situation is more involved since one can allow for 
large values of Aq with some fine-tuning of the uiq/Aq ratio. The relations 



corresponding to pi) and (35) are respectively 



- 1 1 / J_ 

9 < 2 (4p - \) V + 6« 

> (V^T^T (36) 



yff^Tf? 2 V 47 (1 + a) 

For moderate values of a, ( |36|) corroborate our conclusions from (|34l ) and 



(35). In the fine-tuning limit, a — ► 0, the singular behaviour of n in the 
parameter a is apparent. However, in this case our approximations are not 
powerful enough to explain the lower bound Q s > 6\/hf + h\ that is found 
in our numerical analysis. 

Another interesting constraint follows from (|20| ) which implies — (m 2 , + 
(mf. + + Af) < /i 2 + M§/2. Neglecting the top-stop radiative cor- 
rections and using (^) and (^|) yields, 

/t 2 = AV > (4p - 1) M 2 (37) 
From (|l7|) and (p6|), one derives the limits: 

(A 2 /k 2 ) > (4p-l) (A /M >0) 

> (4p - 1) (3a/(l + 2a)) (A /M < 0) 1 ; 

For moderate values of a this implies k 2 < A 2 . At the price of some fine- 
tuning and Aq/Mq < one can obtain k slightly larger than A. 
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Concerning the important parameter tan (3, some simple qualitative pre- 
dictions can be obtained from the following expression, 



^ , = sin 2/3 (39) 

m\ + m\ + ( + + A\ + fi 2 J [5 t L + 2p? 

which is equivalent to ( p0[ ) by the minimum conditions ([H])- Notice that 
sign(tan (3) = — sign(Z?/t). From the renormalized expression (55) in the 
appendix, one sees that for AqMq < the sign of B is preserved and its 
magnitude increased by the gauge renormalisation proportional to Mq which 



increases the numerator in (3S). Therefore, in this case, tan/3 > and it can 



take relatively small values for \Aq\ > \Mq\ . Instead, for AqMq > 0, where 



(27) also applies, the Mq term in (|55| ) tends to exceed the Aq term so that 
in most cases tan (3 < and its magnitude is relatively large. (Exceptions 
are the special cases with small p and |^o| > 2|Mq| with a positive large 
tan fi.) 

Therefore, with these simple analytic approximations for the vacuum 
conditions, together with the approximated condition from the experimental 
limits on the gauginos (to be discussed later on), one easily understands the 
qualitative pattern of the allowed parameter space as obtained from the 
detailed numerical investigation. For instance, (|35"|) and (|36|) are related 
to the decoupling of the singlet sector, whereas sizeable mixing - which 
requires s/ \/h\ + h?, ~ 0(1) - is possible only with the fine-tuning (and 
signs) discussed below d36|). 



6 Qualitative aspects of the mass spectrum 

Let us first consider the chargino masses given by the expression 



(40) 



± \ J [Ml +p? + M*rY- (2£M 2 - sin 2/3) ' 



with M^/Mq ~ 2/3. For the LEP experimental limit we take rn x + > 65 
GeV. Of course, the constraints from this limit on the parameters are the 
same as in the MSSM (with tan /3 in the more restricted range determined by 
the (M+1)SSM as discussed in section ||). On the other hand, this constraint 
plays an important role in the discussion of the previous section as well as 
in this section. 

Because the scale of the soft supersymmetry breaking parameters is set 
by the weak scale, we conveniently express the dimensionful parameters in 
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terms of Mz in this section. One can use the experimental bound, together 
with (p0|) and the (|3ll), to deduce the relations 

A 2 > Up -f\Mi-^> Mf > m 2 x+ > M 2 z /2. (41) 
Let us now define a parameter a by, 

A 2 + ^ = (a-l)M| (42) 

which, for tan f3 > 2 (which turns out to be always the case in our numerical 
analysis), satisfies the approximate inequalities 6/9+1/4 < a < 6p+2. Then, 
from © and © one gets the following expression for the lightest chargino 
mass, 



2 (a- 1) ,, 2 A 3M 
i 2 + ^ '-Ml; - ( 1 + ] 

x a 

so that the approximate bounds 



m J + = ^M|- I 1 + ^1^, (43) 



— M s 2 --M| ~ -M?--Mj- < m 2 + < -Ml --Ml ~ — M, 2 --M| (44) 
14 3 2 z 4 2 3^~x+~8 2 6^ 25 3 6 

are obtained by varying the parameters p and tan/?, and by using the re- 
lation between the gluino and wino masses M|/M| — 10.7. This defines a 
relatively narrow band for m x + as a function of the gluino mass M3, that is 
reasonably reproduced in our numerical analysis. It just expresses the fact 
that the lightest chargino is mostly a wino, with a component of 0(1/ a) of 
the higgsino. 

Now, from ([ij]) and the experimental LEP bound, ""t- 2 + > M^/2, one 
gets a limit on the gluino mass, M3 > 3Mz = 273 GeV (from the previous 
LEP bound, m x + > Mz/2, the corresponding result was M3 > 180 GeV, 
but needed a more detailed approximation than (fi"4])). This corresponds 
to Mq > 100 GeV. Notice that these are only approximations aimed to 
explain rather than predict the numerical results that take all the parameters 
into account. As already stressed, this result is analogous to the limits on 
the chargino masses in the MSSM. For Mq ^> Mz the off-diagonal terms 
become relatively unimportant and fx 2 — » (a — 1)M|. Therefore the lightest 
chargino is mostly a wino with m + ~ M2 and the other chargino has mass 

Xi 

m+ ~ \fa — lm + > \^m+. 

X-2 Xl Ai 

We next consider the neutralino sector, with five Majorana fermions. It 
is particularly important since it includes the lightest odd R-parity state, or 
LSP, and since its phenomenology can deviate from that predicted in the 
MSSM because of the presence of the singlet fermion (singlino, s). The mass 
terms for the neutralinos are, 
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2vs + =M Z ( hi sin \ + /?o cos ; 



+ jxhih 2 + M\BB + M2W3W3 
+M Z (b sin d w + W3 cos O (hi cos /9 + h 2 sin 0\ (45) 



The point to be stressed is the relatively small mixing of the singlino to 
the higgsinos. This is because s 3> /12 > h\ and (2/3) |^4o| < |12i/| < |^4o I > 
fi = As > \/2-^o > y/2Mz- Hence the singlino remains an almost pure state 
of mass \2u\, and the higgsino-gaugino sector is analogous to the charged 
one. The experimental constraints on the chargino imply, as just discussed, 
that the lightest chargino is mostly a gaugino, which will also imply an 
analogous situation in the non-singlet neutralino sector. 

One finds indeed that the lightest states are ^ one which is mostly B, 
with a mass ~ M\ ~ Mo/\/6, and another one, mostly W3, with mass ~ 
M2 — 2M%. Prom the bound on the chargino mass one has M\ > M|/5. The 
most important difference between the (M+1)SSM and the MSSM could be 
in the LSP character: unless one badly violates the universality assumption 
for the scalar masses, the LSP is mostly a B in the MSSM, while the s is 
also a LSP candidate in the (M+1)SSM. Let us consider the condition for a 
singlino LSP, 

\2is\ < Mx (46) 
which in terms of the bare parameters yields: 

\A \ < M ^ < */f Mo (47) 

3 + VI + 8a V 8 y ' 

This corresponds to some constraints in the coupling constants as well. 
From (|34]) one infers n/g < .01. The S v.e.v. has to be relatively large 

because of (^), which gives s/ L/K{ + k£) > 20. From © and © 

one 

derives the constraint, \lijv\ = X/k > 5. On the other hand, (47) also implies 
Mq > 24 ?71q. In other words, the light singlino scenario corresponds to 
gaugino masses being the largest soft terms, a "gaugino dominated scenario" . 
This implies many correlations between sparticle masses. 

If the singlino is the LSP, the second lightest neutralino being mostly a B 
should be quite long-lived. Indeed, its decay into s has to go through its mix- 
ing to higgsinos. Then the decay B — > J is proportional to g\ {h\ + /i^) / s 2 < 
gf/400. Conversely, if (|47|) is not realized, then s will decay into the LSP, 



3 As a matter of fact, when the mass spectrum of the model is relatively low, the mixings 
are more important in the neutralino sector, and, besides the singlino, the lightest state 
can be more like a photino (7) than a "bino" (B). But this can be easily incorporated in 
the discussion that follows where a B-state is assumed for simplicity. 
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mostly B, with a similar coupling (s would be produced in Higgs or higgsino 
decays) . 

Let us now turn to the slepton sector. The lowest lying states are the v 
sneutrinos V e and the "right-handed" sleptons, £r, I = e, /i, r 0. This follows 
from the assumption of universality for the soft terms. The approximated 
expressions for the slepton masses as given in the Appendix are 



m~ = ml + .15M 2 + M§ sin 2 8 w \cos2/3\ 



ml = m z + .5M$ - |cos 2/3| (48) 



ml + .5M 2 - -M% |cos2/3| 



Searches at LEP 1 put a limit of roughly Mg/4 on both of them. For 
the relatively large values of tan/3 that tend to prevail in the (M+1)SSM 
I cos 2/?| ~ 1. With the LEP 1 results on the chargino masses, as pointed 
out above, one had a bound 2Mq > M| and the bounds on m~ gave some 
useful information. With the LEP 1.5 results, the limit on Mq implied by 



(45) is enough to ensure > mj- . 

A more interesting issue is provided by the interplay between chargino 
and slepton searches at LEP 200. In order to compare their masses, we first 
replace Mq in (HI) in terms of the lightest chargino mass in (|43|) to obtain 



2 a 2 2 Ml / 3 3 \ . 2 2, „ /in . 

m k * i^TTj (1 + Ta + a ™S + 7 MI (49) 

which gives a qualitative explanation of figure (2) in the next section. The 
difference between m^r and m x +/2 is almost a measure of the parameter 
mo in the future searches for these particles at relatively high energies. As 
a matter of fact, this reverses the situation at LEP 1, where the chargino 
was expected to be lighter than the selectron. Indeed by comparing ([48]) 
and fl43|), we find that the chargino is always lighter than the selectron if 



M 2 2 < Mf/2, or, equivalents, as m 2 + < M§/3 = (52.5GeV) 2 . This is quite 
well reproduced by our numerical results. 

In particular, in the light singlino scenario, the term ml can be almost 
neglected in (|^) so that the selectron and chargino masses are strongly 
correlated. 



4 In a GUT context, the r superpartners can become considerably splitted and much 
lighter than the other sleptons, as a consequence of renormalization between the Planck 
and GUT scales, as noticed in pi| . However in this paper we are assuming ht > h T at 
Mpi anc k and we consistently neglect this possibility. 
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7 Numerical analysis, experimental and theoreti- 
cal constraints 



In this section we describe the numerical procedure employed to produce the 
figures (1) - (6). As mentioned above, the essential independent parameters 
of the model are (apart from the lepton and light quark Yukawa couplings 
and the known gauge couplings) three Yukawa couplings ht, A and k and, 
under the universality hypothesis at Ao = Mqut, the three soft parameters 
Mo, mo and Aq. Since these are the only dimensionful parameters in the 
theory and the scale is finally set by the experimental value of Mz, the true 
independent parameters are the five dimensionless quantities h to , Ao, Ka, 
m /M and A /M . 

The numerical results presented here have been obtained by scanning 
over ~ 10 6 points in this five dimensional parameter space. (Here we used a 
logarithmic measure, but we have checked the results with different measures 
as well. A lot of scanning of particular regions of the parameter space has 
been carried out in order to verify the boundaries of the parameter space 
and the mass ranges, but the results of these scannings are not shown in 
the figures.) In each case we integrate the one loop renormalization group 
equations down to the electroweak scale v = 174 GeV. To one loop accuracy 
it is sufficient to use v as the low energy scale, independently from the exact 
particle masses of 0(y). At the one loop level, one can also take advantage 
of (H) and compute the coefficients in these equations, what saves a lot of 
computer time to be used for a better scanning. 

Actually, also two loop renormalization group equations have been con- 
sidered in [Q]. In this case, however, the low energy scale has to be defined 
more precisely: The decoupling of particles with masses of 0(v) has to be 
properly taken into account; these masses, in turn, are only known once 
the Higgs potential (which depends on these masses via the radiative cor- 
rections) has been minimized. As a consequence, the numerical procedure 
becomes much more envolved and allows only to study much less points 
in the parameter space. On the other hand, the corresponding numerical 
results hardly deviate from ours. Individual points in the parameter space 
may well lead to somewhat different particle masses, once two loop contribu- 
tions have been taken into account. Allowed mass ranges and correlations, 
however, remain practically unchanged. We found it much more important 
numerically to minimize the full one loop Coleman- Weinberg Higgs poten- 
tial including the non-logarithmic contributions than to include two loop 
logarithms. 

Having obtained the parameters at low energy it is most convenient to 
check first the absence of slepton v.e.v.'s according to eq. ([24]) at the scale 
~ A e /h e . Next we minimize in each remaining case the effective potential 
(H) including the Coleman- Weinberg radiative corrections (H|),(10), numer- 
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ically. We test, whether the minimum with (S), (Hi) and (H2) 7^ is the 
lowest one, and dismiss the corresponding set of parameters otherwise. We 
calculate tan (3 and discard the initial conditions that lead to tan (3 > 30 
(since we assume ht 3> /if, in the present work; larger values turn out to be 
somewhat disfavoured in the model). 

In the remaining cases we determine the overall scale of the dimensionful 
parameters by identifying (Hf} + (-fff ) with 2M^/(gf + g^), and compute 
the physical masses of all particles. Then we impose the following exper- 
imental constraints: For the top quark pole mass mt (where we take the 
leading QCD corrections to the pole mass into account |28|]) we require 
168 GeV < mt < 192 GeV |J. Next we impose the LEP 1.5 lower bound on 
the lightest chargino mass m x +: m x + > 65 GeV [jTj. As a result, nearly all 
other experimental bounds on new particles turn out to be satisfied auto- 
matically, with some exceptions discussed in section || Indeed, though the 
experimental bounds on sleptons, stop, gluinos, as well as on the decays of 
the Z into neutralinos and Higgs scalars are implemented in our code, these 
conditions only eliminate a few marginal points in the parameter space. This 
remarkable property also can be seen from the plots in the next section. 

Let us briefly discuss the range of the bare parameters, which turns out to 
be consistent with our theoretical and experimental constraints. Concerning 
the Yukawa couplings we first remark that, as in the MSSM, the known range 
of the top quark mass together with the solution of the RGE for ht as given 
in the appendix allows easily to obtain the allowed range for ht : mt > 168 
GeV leads to ht > 0.468. The other two bare Yukawa couplings Ao and kq 
turn out to be fairly small, 2.7- 10~ 3 < A < 0.32 and 1.1 • 10" 4 < k < 0.33, 
and relatively closely related, 0.04 < kq/\q < 1.14. 

As already remarked below eq. (0), the experimental lower limit on 
m x + implies a lower limit on the soft susy breaking parameter Mq. The 
numerical analysis is in good agreement with the approximate analytic result 
and leads to Mq > 90 GeV. No strict upper limits (in the absence of fine 
tuning constraints) on the soft susy breaking parameters have been obtained. 
The scalar mass term m\ turned out to be unconstrained by our analysis; we 
assumed, however, m\ > 0. The numerical analysis also confirms the allowed 
range for the third susy breaking parameter Aq as a function of Mq and 
mo, eq. (28). The corresponding lower limit on Aq is a particularity of the 
(M+1)SSM, since it is required by the need to destabilize the scalar potential 
in the singlet direction. The correlations between the particle masses and 
the bare parameters can finally be obtained from the approximate analytic 
solutions given in sect (6) and the appendix; in table 1 we present the precise 
numerical parameters and masses for two particular generic cases: one with 
Aq negative and a non-singlet LSP, and one with Aq positive and a singlet 
LSP. 
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8 Predictions for the particle spetrum 



In order to allow to study the correlations among the particle masses we 
choose to plot all masses against the lightest chargino mass ni x +. Each 
point in these plots corresponds to one of the ~ 1.5 ■ 10 3 points in the 
parameter space, which satisfy all our theoretical and phenomenological 
constraints. The density of points in the plots is clearly not uniform; regions 
of low density (if not completely empty and hence forbidden), typically 
towards larger particle masses, correspond to regions where more and more 
fine tuning is required. We will not, however, employ quantitative "fine 
tuning constraints" in this paper in order to constrain the particle masses 
from above; a reader, who is interested in "probable particle masses", can 
obtain them by investigating the relative densities in our plots. Let us first 
consider the gluino mass M3, fig. 1. As already noted in eq. ([44|), given 
the present experimental bounds, m x + and M3 turn out to be strongly 
correlated. Fig. 1 agrees well with the approximate analytic relation (|44]), 
and one sees that the lower experimental bound on m x + turns into a lower 
bound M3 > 260 GeV on the gluino mass within our model. For large 
masses, the ratio between the lightest chargino and the gluino is close to 
92/93 as P r edited by the chargino being a wino, but slightly smaller as the 
parameters in the chargino mass matrices are dynamically correlated (as 
discussed in section |6|). 

Next we turn to the slepton masses. The lightest charged "right-handed" 
slepton mass mj- is plotted against m x + in fig. 2. We see that m x + > 65 
GeV implies mj > 55 GeV in our model. In most cases, as noted below 
(p9|), the charged slepton mass m j R has the tendency to be smaller than 
m x + (except for large values of the bare scalar mass mo), although it always 
satisfies fnj R ^ m x + m agreement with the analytic approximation (|4^) . 
In fig. 3 we plot the sneutrino mass my against rn x + • Given m x + > 65 GeV, 
the LEP 1 lower limit on only eliminates a tiny additional corner in the 
parameter space. 

The neutral scalar Higgs sector has already been discussed in detail in 
ref. ||, in particular with respect to future experiments at LEP 2. Although 
in m only the LEP 1 bound on the chargino mass rn x + had been taken into 
account, the results concerning the neutral Higgs sector remain essentially 
unchanged. Therefore we briefly repeat only the essential features. In order 
to take the experimental constraints on neutral Higgs bosons properly into 
account both their masses and their couplings to the Z bosons have to be 
evaluated. In fact, the negative LEP 1 results on neutral CP-even Higgs 
scalars || eliminate only a tiny region of the parameter space of the present 
model, corresponding to rrih < 58 GeV and a coupling hZZ as large as in the 
non-supersymmetric standard model. On the other hand, a large region with 
smaller values of m^, but small couplings hZZ remains unconstrained by 
LEP 1. In this region the lightest CP-even Higgs scalar is essentially a gauge 
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singlet and hence decoupled. In view of this possibility the experimental 
searches for and the theoretical upper limits on the masses of the lightest 
CP-even Higgs scalar have to be turned into considerations of the "lightest 
visible CP-even Higgs scalar". Fortunately it turns out within the present 
model |, |, ||, that upper limits exist also on the lightest visible CP- 
even Higgs scalar (which could be the second lightest), varying from 140 to 
160 GeV for gluino masses (as a measure of the susy breaking scale) from 
1 to 3 TeV. Similar complications arise also in the CP-odd Higgs sector, 
which contains a gauge singlet state as well. Here no upper limit can be 
derived, and we just remark that a visible CP-odd Higgs scalar below 130 
GeV is impossible within the present model and the present experimental 
constraints on m x +. 

As in the MSSM the mass of the visible (non-singlet) CP-odd Higgs 
scalar is actually just somewhat below the mass mjj+ of the charged Higgs; 
cf. the upper left entry of the CP-odd mass matrix eq. (15) and eq. (16). 
In fig. 4 we plot + against rn x + , and we see that the present lower limit 
on m x + implies m H + > 160 GeV (increasing rapidly with increasing lower 
limits on m x +). 

The range of tan (3 is fairly restricted, and depends strongly on the sign 
of the bare susy breaking parameter Aq: for Aq > we have tan/3 < —2.6, 
whereas for Aq < we find both signs, but large absolute values for tan/3: 
tan/3 < —6.5 or tan/3 > 8.7. Since these inequalities hardly depend on the 
mass of the chargino, a figure does not provide additional information. 

Next we turn to the squarks. Because of renormalization effects between 
Mqut and Mz (cf. the appendix) the lightest stop eigenstate is the lightest 
of all squarks; subsequently we concentrate on this particle. In fig. 5 we 
plot its mass against rn x +, and we see that the present lower limit on 
m x + implies ma\ > 170 GeV. A possible "supersymmetric" top quark decay 
mode is thus practically excluded within the present model. 

Finally we consider the neutralino sector. As discussed in chapter (6), 
two different scenarios are possible within the present model, depending on 
the mass of the nearly pure singlino: In the case of a heavy singlino, the two 
lightest neutralino states are a practically pure bino with a mass Mi and 
a practically pure wino with a mass M2 ~ 2M\. The numerical procedure 
confirms this approximate analytic result to a very high accuracy. Since 
furthermore the chargino mass rn x + is close to M2, cf. eq. (f44|) , the two 
lightest neutralino masses can easily be obtained as a function of rn x + and 
a figure is of no further use. Here the LSP is, of course, the bino, and the 
supersymmetric decays of all sparticles will proceed as in the MSSM with a 
suitable set of parameters. 

The case of a light singlino is actually more interesting. From the numer- 
ical analysis we find, that now the second and third neutralinos correspond 
to the above mentioned nearly pure bino and wino states. In fig. 6 we plot 
the mass of the lightest singlino, m/ s , against the bino mass M\ (for those 
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parameters where mi s < Mi). Now supersymmetric decays of all sparticles 
will start as in the MSSM, until the bino state is produced. (A direct decay 
into the singlino LSP will be heavily suppressed due to its tiny coupling.) 
Finally the bino will decay into the singlino LSP, producing an additional 
cascade, whose consequences on experimental susy searches still have to be 
worked out in more detail. 



9 Conclusions and outlook 

In the present paper we have presented an exhaustive overview over the 
parameter space and the particle spectra, which are consistent with the 
correct vacuum structure and present experimental constraints, within the 
(M+1)SSM under the universality hypothesis. Since the scalar potential 
and the corresponding minimization conditions differ considerably from the 
MSSM due to the presence of the additional singlet superfield, the outcome 
of this analysis was far from predictable. We have seen that still a large 
region in the parameter space consistent with all constraints exists. This 
region is characterized, however, by a) small values of the new Yukawa cou- 
plings A and k, b) a large value of the scalar singlet vev, c) a strong tendency 
of the singlet states to decouple from the rest of the spectrum. In spite of 
the complexity of the model this region allows for fairly accurate analytical 
approximations to the vacuum minimization conditions and particle masses, 
which we have worked out and presented in some detail. 

The near decoupling of the singlet states leads to a particle spectrum, 
which is close to the MSSM for a suitable set of the corresponding parame- 
ters. Our figures show strong correlations or inequalities among the particle 
masses, which allow for experimental falsifications of our underlying hy- 
pothesises. For instance, as in the MSSM with the universality hypothesis, 
sparticles are not light enough to play an important role for the radiative 
corrections to R b i [Q. Only the Higgs boson, the lightest chargino, sleptons 
and neutralinos could be light enough to be discovered at LEP 2; gluinos, 
squarks and charged or CP-odd Higgs scalars are necessarily too heavy in 
our model. 

Our approximate analytic results allow to translate future possibly nega- 
tive sparticle searches into constraints on the soft susy breaking parameters. 
Obvious are lower limits on the universal gaugino mass Mo from negative 



chargino searches (cf. eq. (44)) and, of course, negative gluino searches. 
From eqs. (50), one can furthermore easily deduce excluded regions in the 
trq — Mq plane from negative "right-handed" slepton and sneutrino searches. 
These relations are actually not different from the MSSM. Eq. (28), how- 
ever, allows to turn these constraints on mo and Mq into constraints on Aq, 
which is a particular feature of the (M+1)SSM. 
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If future sparticle searches turn out to be successful, and the mass pat- 
tern predicted by our model turns out to be confirmed, the question arises 
whether the model can be distinguished from the MSSM with correspond- 
ing parameters. The answer depends on the corresponding parameters of 
the (M+1)SSM. First, there still exist small corners in the parameter space, 
where mixing with the singlet states either in the neutral Higgs sector or 
in the neutralino sector could be possible; in this case the couplings of the 
corresponding visible particles could be substantially and hence measurably 
smaller than in the MSSM. A more likely possibility - but by no means cer- 
tain - is the singlino LSP scenario. As discussed in the previous section, in 
this case sparticle decays will produce an additional cascade with respect to 
the MSSM. This scenario desserves further phenomenological investigations. 
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A Appendix 



In this Appendix we list the analytic solutions of the RGE for soft susy 
breaking terms in the approximation where the dependence on all Yukawa 
couplings but the top one, /if, are neglected. 

Also, flavour mixing is neglected in these expressions (see [^]) for solu- 
tions that fully take family mixing into account). However, for the sake of 
generality, universality is not assumed for the soft terms at the unification 
scale, Ao- All quantities are defined at the scale A (unless stated other- 
wise) and the boundary conditions at the scale Ao, are generally denoted as 
Xq = X(Aq). We define the following quantities: 



1 



167T 2 A 

,2 



no, = U(f) 

a= i \ya / 

t 



J = ^(f)J U {u) (t')dt' (50) 

The index a = 1,2,3, refers to the U(l), SU(2) and SU(3) gauge groups, 
respectively, and Cq, is the sum of the their quadratic Casimir eigenvalues 
for the three fields in the i-th Yukawa coupling. We replace the integral J 
by the value h CT n of ht at the scale A corresponding to a Landau pole at Ao, 
by defining 



J- = 12 J (51) 

' l crit ' l crit 

The solutions for the Yukawa coupling are then as follows 



A 2 


= A§n (A) (i- 


-pi 12 


hi 


= ^ 2 n H (i 


-p) 


2 

■c(u) 


- h 2 t ^ n (u) 


(i-p) 1/2 


h 2 b 


= h 2 b0 U (d) (l 


-p) 1/6 


.2 

's(d) 


= h2 s(d) U (d) 




hf 


= KP-M 


{£ = e,fi,r 


K 2 


- K 2 





(52) 

The soft trilinear terms are as follows: 
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Ax = A Xo --p(A to +ZM ) + 2M t(3g% + gj) 
A K = A k0 

(16 13 \ 

1 16 13 \ 

= A c{u)0 - -p (A t0 + £M ) + 2M t ( —gl + 3# 2 2 + y y 2 J 

A b = A m -^p{A to +tM Q ) + 2M Q t(^g 2 + 3g 2 + ^g 2 ^j 

A s{d) = A s{d)0 + 2M t(^-g 2 + 3g 2 + p 2 ^ 
A tM = A T ^ e)0 + 2M t(3g 2 + 3g 2 ) 

Z = (j- 1 ) = (I2^ ri t* - 1) (53) 

Soft supersymmetry breaking scalar masses can also be written in terms 
of p, t and J (or £), as follows: 

m 2 s = m| 

m \ = mj + gftso + 2j L tMQ 

m \ = ^20 - y ™o ~\ K ~ 9i ts o + 2-ihtMl 

m\ = mf- - pml - + ^ 2 is + 2-y v tMQ 

m Q3 = m Q30 ~ |"*o - I K - gflfrao + 2^ Q tMl 

m 2 B = m 2 B0 - ^gfts + 2^ D tM^ 

m u t = m Uio + ^9i ts o + 2~iutMl 

m Qi = m Qio ~ + 2-iQtMl 

m D t = m Dio ~ g0i**o + 2~i D tMl 

m Lj = m LjO + 51*50 + 2j L tMQ 

m Ej = m Ejo - 2glts + 2~/ E tMQ 
7 l = 3g% (l - g%t) + <fi (l - ll(fit) 

1U = yff 2 2 (l + 3 5 2 t) +^ 5l 2 (l-ll<7 2 t) 

16 1 
7Q = -gl{l + 3g 2 t)+3g 2 {l-g 2 t) + -g 2 {l-llg 2 t) 
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ID 
IE 

SO 



ml 



K 



y gl (i + 3 5 f t) + \g\ (l - n^t) 

4 ffl 2 (1 - ll 5l 2 t) 

2 2 i / 2 , 2 , 2 2 o 2 \ 

™ff 2 - niHi + 2^ v m Qi + m A + m Ei ~ m L l ~ 2 ™Ui) 

i 

1 / 2 , 2 , 2 \ 
g { m 2 + m T + m Q) t=0 

(i-p) (^ + eM ) 2 -^M 2 + 

+2 (y fli + 3 52 2 + y 5l 2 ) (^ + l)tM 2 



t=o 



(54) 



Finally, let us recall that the renormalization effects on the effective 
variables B and fx are easily obtained by taking the limit in which the singlet 
fields are decoupled at its classical value. Then, 



B 



A x + is = B --p (A t0 + £M ) + 2M t (Sgj + gf) 



p = aV = AgnW (i - P f 2 



(55) 
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Figure captions 



Figure 


1: 


Gluino mass M3 versus the chargino mass m x + (all in GeV) 


Figure 


2: 


Slepton mass rnj R versus m x + 


Figure 


3: 


Sneutrino mass versus m x + 


Figure 


4: 


Lightest stop mass versus m x + 


Figure 


5: 


Charged Higgs mass m H + versus m x + 


Figure 


6: 


Mass of the second lightest neutralino m x o versus the mass of 



the lightest neutralino m o in those cases where xi 1S essentially a 



singlino 

Table caption 

Table 1: Parameters and particle masses (in GeV) for two particular cases 
(with M susy not too large). The three neutral CP even Higgs scalars 
are denoted by h\ , /12 and /13 (in the order of increasing masses) , the 
two neutral CP odd Higgs scalars by p\ and P2, and the two lightest 
neutralinos by x? and respectively. In the first case x? is essentially 
a non-singlet, in the second case it practically a pure singlino. 
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